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Abstract 

The three-body continuum Coulomb problem is treated in terms of the generalized parabolic 
coordinates. Approximate solutions are expressed in the form of a Lippmann-Schwinger type 
equation, where the Green's function includes the leading term of the kinetic energy and the total 
potential energy, whereas the potential contains the non-orthogonal part of the kinetic energy 
operator. As a test of this approach, the integral equation for the (e~, e~, He"^'*') system is solved 
numerically by using the parabolic Sturmian basis representation of the (approximate) potential. 
Convergence of the expansion coefficients of the solution is obtained as the basis set used to describe 
the potential is enlarged. 
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I. INTRODUCTION 




The three-body continuum Coulomb problem is one of the fundamental unresolved prob- 
lems of theoretical physics. In atomic physics, a prototype example is a two-electron con- 
tinuum which arises as a final state in electron-impact ionization and double photoion- 
ization of atomic systems. Several discrete-basis-set methods for the calculation of^uch 
processes have recently been developed including convergent close coupling (CCC) 
the Coulomb-Sturmian separable expansion method 4], the J-matrix method pl 
all these approaches (see also [sl) the continuous Hamiltonian spectrum is represented in 
the context of complete square integrable bases. Despite the enormous progress made so far 
in discretization and subsequent numerical solutions of three-body differential and integral 
equations of the Coulomb scattering theory, a number of related mathematical problems 
remain open. Actually, the use of a product of two fixed charge Coulomb waves for the 
two outgoing electrons as an approximation to the three-body continuum state is typical 
of these approaches. As a consequence, a long-range potential appears in the kernel of the 
corresponding Lippmann-Schwinger equation. Since this integral equation is non-compact, 
its solution therefore is divergent as the size of a basis used to describe the potential is 
increased. Note that a renormalization approach can not cure this problem. 

A theoretical treatment of Coulomb breakup problems (which does not require screening 
or any regularization) has been recently suggested by Kadyrov et al. [9]. 

In several papers 10l-ll2| a new approach for solution of the three-body continuum 
Coulomb problem was introduced. The development of the method is chiefly based upon 
the fact that the asymptotic wave operator, which determines the wave function behavior 
when all interparticle distances are large (in the Qq domain), is separable in terms of gener- 
alized parabolic coordinates 13(|. The parabolic coordinate eigenfunctions of the asymptotic 
wave operator which satisfy the Redmond's conditions (l4] in the asymptotic domain 
are expressed in terms of a product of three two-body Coulomb functions, each depending 
of one parabolic coordinate. These functions, often called the C3 (or BBK) wave functions 



ISMIT]. are successfully used as the final-state wave functions for calculating cross sections 
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20|. In it has 



for electron-impact ionization and double photoionization of helium 
been proposed to use an integral equation of a Lippmann-Schwinger type to construct an 
approximate solution that describe three charged particles moving in the three-body con- 
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tinuum. In this integral equation the Green's function in the whole configuration space 
(and not just in the domain Qq) is approximated by the inverse of the asymptotic wave 
operator. In turn, the well-known non-orthogonal part of the kinetic energy operator, which 
represents the difference between the total and asymptotic wave operators, plays the role of 
the potential. Asymptotic behavior of solutions is determined by the inhomogeneous term 
which is given by the C3 wave function. 

To test the practicality of the Lippmann-Schwinger equation approach, we restrict our- 
selves to the so-called outgoing approximation [2l|, which assumes that the sought-for solu- 
tion (as well as the inhomogeneity) depends only on the parabolic coordinates j = 1, 2, 3, 
and therefore the approximate potential operator contains only the terms which involve the 
mixed derivatives d"^ /dC,idC,s, I ^ s. In this work we address the question of compactness of 
the kernel of the integral equation. Any compact operator may be approximated arbitrary 



closely by an operator of finite ran! 



^. In order to study the properties of the kernel we choose 



the parabolic Sturmian basis set 22| and construct a sequence of separable kernels. Then 



to check the existence of a limit of this sequence we examine the convergence behavior of 
the first few expansion coefficients of the solution as the basis set is increased. The coef- 
ficients are found to exhibit oscillations whose amplitude does not decrease as the number 



of terms in the representation of the potential grows. T 



lis result is similar to the Gibbs' 



23|, 



24|). where the oscillation of an 



phenomenon known from Fourier analysis (see, e. g., 
approximant about the exact function (which possesses a discontinuity) is a consequence of 
the abrupt truncation of the Fourier sum. In order to avoid or at least to reduce the Gibbs' 
phenomenon, smoothing procedures are used that attenuate the higher order coefficients 



23l . |2J]. In this paper, we use the Lanczos smoothing factors, introduced in the potential 
separable expansion (PSE) method \^ (see also 26| and references therein), in constructing 
the basis set representations of the potential. 

This paper is organized as follows. In Sec. II we introduce the notations, recall the gen- 
eralized parabolic coordinates definition and express a formal solution for the three-body 
Coulomb problem in the form of the Lippmann-Schwinger-type equation. In Sec. Ill we 
briefiy outline the parabolic Sturmians approach. In particular, we present the matrix rep- 
resentation of the three-body Coulomb Green function and consider the potential operator 
approximation. In Sec. IV calculations of the continuum state of the (e~, e~. He''"'') system, 
where both electrons recede from the residual ion in opposite directions with equal energies, 
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are described. Our aim is to study the rate of convergence as the basis set used to describe 
the potential operator is enlarged. If the same number of basis functions for each parabolic 
coordinate is used, the problem of numerical solution rapidly gets out of hand. Thus, in 
the separable expansion of the potential, the number of basis functions for the three chosen 
curvilinear parabolic coordinates ^i, ^2, ^3 is increased, whereas a single basis function is 
taken for each of the remaining three coordinates rji, 772, V3- The calculations show that the 
convergence on a basis of reasonable size can be obtained by using the Lanczos smoothing 
factors. Sec. V contains a brief discussion of the overall results. Atomic units are used 
throughout. 



II. THEORY 

Wc consider three particles of masses mi, mi, ms, charges Zi, Z2, Z3 and momenta ki, 
k2, The Hamiltonian of the system in the center of mass frame is given by 

1 . 1 . Z-i Zo Zi Z-x Z-\ z^ , , 

E = -^An - — Ar + + + 1 

2//12 2/^3 ri2 r23 ri3 

where ris denotes the relative coordinates 

ris = r/ - r^, = \ris\ , (2) 
R and r are the Jacobi coordinates 

miri + mara , . 

R=ri-r2, r = r3 ■ . (3) 

mi + m2 

The reduced masses are defined as 

mim2 ms (mi + m2) 

/^12 = , 1^3 = ■ 4 

mi + m2 mi + m2 + ms 

In the Schrodinger equation 

= (5) 

the eigenenergy > is given by 

2/^12 2/^3' ^ ^ 

where K and k are the momenta conjugate to the variables R and r. Substituting 

$ = gj(K.R+k.r)^ 



into ([5]), we arrive at the equation for the reduced wave function \I/ 
1 



1. i „ i -. „ ZtZo ZoZ-i ZiZ'i 
2a^12 2/^3 /ii2 /is ri2 r23 ria 



^ = 0. 



(8) 



Leading-order asymptotic terms of ^ in the Qq domain are expressed in terms of the gener- 
ahzed parabohc coordinates 



^1 — 1^23 + ■ r23, rji — r23 — k23 ■ r23, 
6 = '^is + ki3 • ri3, r]2 = ri3 - ki3 ■ rig, 
6 = n2 + ki2 ■ ri2, ?73 = ri2 - ki2 ■ r^, 



(9) 



where k/g = is the relative momentum, k/s = kis = |kjj_^ The operator in the 



square brackets, denoted by D, can be decomposed into two terms 13 1 



D = Do + D,, 



(10) 



where the operator Dq contains the leading term of the kinetic energy and the total potential 
energy: 



^0 = E 



for j I, s and / < s, 
f d d d 

-2(— —-ik —1 



h 



Vj 



(11) 

(12) 
(13) 



Here tis = ^ 



mi+m. 



. The operator Di represents the remaining part of the kinetic 



energ y | l3| which in the case of the (e , e , He"*""*") = (123) system with ms = oo takes the 
form [21] 



+ u7 ■ u 



02 



3 drjj 8^3 



+ u7 • u. 



+ d2 

■ ^^3 a^, dm 



92 



3 drij dr]z 



where 



(14) 



(15) 



The asymptotic behavior of solutions \l/ is determined by the operator Dq. In particular, 
there exist solutions to the equation 



D^^C3 = 0, 
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(16) 



which satisfy the Redmond conditions in Qq. These solutions are well-known the C3 wave 
functions. \E'c3 is expressed in terms of a product of three Coulomb waves. For example, 
\E'c3 with pure outgoing behavior is written as 

3 



^C3 = Y[iFi{ttis, 1; -tkisQ. 



(17) 



n turn, Di is regarded as a perturbation which does not violate the asymptotic conditions 



3- 



Our goal is to construct an approximate solution ^E' of ([8]) that satisfies the boundary 
condition fll7p in the asymptotic domain. For this purpose, we rewrite in terms of 
the operators 



and 



'H= Yl i^j + rjj) Do = /ii3 (^2 + V2) /ii2 (6 + ^3) ^1 

+/i23 (6 + Vl) 1^12 (6 + ^73) h + 1^23 (^1 + Vl) 1^13 (6 + ^72) 



3=1 



after multiplying on the left by Yl f^is {^j + Vj)'- 



n + v ^ = 0. 



(18) 



(19) 



(20) 



(21) 



Thus, given the Green's function operator Q = "H^^, we can take into account the non- 
orthogonal term Di of the kinetic energy operator (which is larger than the total potential 
in the "inner zone" j27|) by putting it into the kernel of the Lippmann-Schwinger type 
equation: 

^ = ^C3 - ^V^. (22) 



Green's functions 

Based on the fact that the original operator Dq is separable in the parabolic coordinates 
([9]), the inverse of the six-dimensional operator "H ( |T8l) can be expressed as a convolution of 
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the three two-dimensional Green's function operators G^^^ whose kernels 



(0, Vj\Gf {tis, £,) % V-) ^ Gf' {Us, £f 0, i'p i) 



(±) 



(23) 



satisfy the equations 

+ ^^ls c, (0 + m)\ of it,., E,- e„ n,, e„ v-) = 5 - e;) 5 {v, - v-) ■ m 

In view of ( ITSl) the separation parameters Cj are subject to the condition 

Ci + C2 + Cs = 0. (25) 



In ( l23l) and we introduced auxiliary variables: 

"2 e 



7l 
2 



— f^ls Cj. 



(26) 



The technique presented in 



28j can be employed to derive useful forms for the two- 



dimensional Green's function. For example, Cj^'' can be expressed in the form (for simplicity 
we omit the indices): 

{t, £■ e, ri, e', V) = e5'=(«'-«+''^''') jdz sinh(;2) [coth (f)]^ 



■^Q±q{i+i'+ri+v') cosh(^) (^i^y^sinh(2)) Jo (Ti7v/Wsinh(z)) , 



(27) 



where 



k 

T=-t, 
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(28) 



/,y(z) is the modified Bessel function of order u |29 1 . 

The resolvent of T-L standard representation 30|] involves integration along two contours 
which encircle the spectra of two of the wave operators f)j, j = 1, 2, 3. In the paper 12| an 
integral representation of suitable for numerical computation has been proposed (see 
below) . 



III. PARABOLIC STURMIANS APPROACH 

If the kernel QV is compact, then the integral equation ( l22l) can be solved by, e. g., the 
algebraic method in which the potential V is approximated by operators of finite rank. For 



22| 



the expansion of V we use a set of square-integrable parabolic Sturmian functions 
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<t>n, m, (0> m) = ^n, (0) , (30) 

ijn {x) = v^e-'^"L„(26,x). (31) 

The basis functions ( 130|) . ( 13T|) are parametrized with a separate Sturmian exponent bj for 
each pair {^j, rjj}, j = 1, 3. Thus, the operator V is represented by its projection V-^ onto 
a subspace of basis functions, 

V^ = ^ ^ 5Z 1^) (^1 V lOT') (^1 , (32) 

■'J •'J 

and the solution \1/ of the problem is obtained for V-^. Inserting V-^ into Eq. fl2^ then leads 
to a finite matrix equation for the expansion coefficients [a\tyi = (^| ), 

a = a^o) - gVa, (33) 

which has the solution 

a= (l + £V)-%(o)_ (34) 

Here [£]f^ryj/ = (^| Q |W) and [V]fyjfy^/ = {^\ V are the Green's function operator and 
potential operator matrices of order A/" = Y['j=i -^j ^^"^ ^^^^ coefficient vector of 

"^cs- The wave function is expressed in terms of the solution of Eq. fl5^ : 

^ = *C3-E E E (35) 

j=l nj=0 mj=0 

where C = V a. 



Green's function matrices 

To construct the six-dimensional Green's function matrix Q_, we need the two-dimensional 
Green's function fl27|) matrix with elements 

oo oo oo oo 

Gi%.,n',m' n = JJJJd^dvdCdr^'{<l^nm\ ^, V) 

(36) 

Inserting (127|) into we obtain after some simple but tedious algebraic manipulations 

M+) (f n = ± f {-0)"+""' „ r(ir+^+i)r(i^+i-2£) 

'~'n,m;n',m'\''^ ) 2'y \ C j f-xW+m r(ir+X+2-£) 

X2F1 {K + l-2e,iT-i; iT + K + 2-i; C~^), 
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where 



min(£, I/) 

Q = 5^ 

j=max(^— /I, 0) 



n\ in 
J J \J 



m 



m 



V = min(n, n), /i = min(m, m'). (39) 



Replacing 7 by —7 {X^9,6^X, (^1/(,t^ — r) in fl37|l gives 



n, m; n' , mJ 



/,m' ^) — ^(C — 1) ^_^^„'+m C£ C 



^=0 



r(-jT+x+2- 



(40) 



X2F1 (is: + 1 - 2£, -IT - -ir + + 2 - £; ()• 



The elements G*'^^ 



Tij , m,j ; n'. , m'j ' 



J = 1, 3 are used in obtaining the matrix Q_ of the six- 
dimensional Green's function operator Q. It has been shown in 12| that, e.g., Q}~^^ can 
be represented in the form of a double integral over the complex variables Si and S2 along 
straight-line paths (see Figure [T]), on which £1, S2 are parametrized by 



where Ei, E2 are real and — tt < (f < 0. Namely, we have 



(41) 



;(+) 



xG 



(27ri)2 /123 Mi3 



00 00 , 

/ / dEidE2Gl';l^.^^,^,[t 



23, ^ + Eie'^ 



—00 —00 



> > (tu; !^ + E2e'n G^^^ , , (tia; ^3). 



(42) 



Here S3 is given by 



'^12 



AH2 



AH2 



^3 = ^- — + — i52 e 



|arg(£^3)| < TT, 



(43) 



,/^23 /^13 

as follows from (l25i) and (l26l) . 

Aside from the replacement V — )■ V^^, we make an approximation, which consists in 
ignoring the correct boundary conditions in two-body asymptotic domains (the asymptotic 
behavior of the C3 wave function in the neighborhoods of the regions Qj, j = 1, 2, 3 has 
been obtained in |31|). 

Further approximations are introduced in the treatment of the potential operator. 
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The potential operator 



As an example of a three-body Coulomb system above the threshold for total break-up, 
we consider a final state for double ionization of helium. Thus, the potential V is given by 
(BU]) and (IT^ . In order to calculate the matrix of V in the basis (12^ . we need to express 
(fT4|) in terms of the parabolic coordinates. From ([9]), it is easy to obtain 

f. ^ _ (e2+^2)' + fe+r/3)'-(ei+^i)' .... 

2 (^1 + Vl) (^3 + V3) ^ ^ 

Whereas, evaluation of matrix elements of scalar products Vij ■ kj^ with {i, j} ^ {/, s] in the 
general case requires the inversion of the transformation (l9l). This (numerical) procedure is 
reduced to finding roots of a quartic polynomial (see, e. g., 211]). Thus, it might appear 
that the corresponding Cartesian coordinates r^^ are complex. To simplify matters, we take 
ki3 = k and k23 = — k and therefore ki2 = k. In this case we have from (|9]) 

ri3-ki2 - — ki3-ri3 - 

?23-ki2 = -;^k23-r23 = -|^, (46) 

ri2 ■ ki3 = — ri2 ■ k23 = " I'l2 = ^3+^3 ■ 

Thus, the matrix of the potential operator V can be constructed in closed form without 
the need for numerical integration. However, it should be noted that in the domains of 
integration fii : + r/i > ,^2 + ^72 + ■Cs + ''?3 and ^2 : 6 + ''?2 > ^1 + ^71 + ^3 + ^73 the triangle 
inequality is violated, so that ri3 ■i\2 < —1 and r23 • ri2 > 1 within these domains. Hence 
the terms ri3 • fi2 and r23 ■ ri2 grow without bound in Qi and ^2, respectively. Obviously, 
such behavior is inconsistent with compactness of the equation. Actually, our calculations 
have shown that the use of f 13 ■ ri2 fj44|) and r23 ■ ri2 f|45l) in the potential operator leads 
to divergence. To avoid this problem, adequate analytical continuations of these scalar 
products into the regions Qi and Q2 should be performed. Here we approximate ri3 ■ ri2 and 
^23 ■ ri2 by their projections in the direction k: 
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Note that at least the absolute values of the terms on the right-hand side of ( H7I) and (HHI) 
as well as of f l46|) . are bounded by 1. Thus we obtain the approximate formula 

+2 (6 + rii) (-6^6^ + 6^ + V2£^^3^ - m£^) ■ 

Further, we take into account only the mixed derivatives d^/d^id^s and 8^/8^29^3- Thus, 
in our calculations, we use the potential 



(49) 



V = -2 

which corresponds to the outgoing approximation 



(50) 



2l|. 



The inhomogeneity 

The expansion coefficients [a'-^-'j^^ of the inhomogeneity \l/c3 (flTl) is written in terms of 
polynomials pn [Uj : 

[^''% = n ^2 {b, + ^k.)r'^ i-irp.^ L + , (51) 

Pnir; = W -n, - + it; -n + - + it; C • (52) 



nl r(i-zr) 
The constant factor 

[^^°^]o = rir[2(&.+^^^^r^'^ (53) 
is omitted below for simplicity, so that we consider the "reduced" coefficients 

[a]^^[aU/[a^X. (54) 



IV. RESULTS FOR THE (e , e , He++) SYSTEM 

Let us consider the case of a back-to-back electron emission with equal energy sharing. 
We put kis = k = 1.5 and choose the values of the exponents bj in the basis to be equal to 
the wave number, i. e., bj = b = 1.5, j = 1, 2, 3. 

We use the single basis function ipQ (!3T|) for the parabolic coordinates rjj and up to sixteen 
functions for each of the three coordinates C,j in the potential operator expansion ( 1321) . Thus, 
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we put Ml = M2 = M3 = 1 and examine the convergence behavior of the first expansion 
coefficient [a]f^ as the number N = Ni = N2 = N3 of the basis functions ipnj (0)? J = 1; 2, 3 
is increased. 

The only practical limitation on the total number A/" = A^^ of the basis functions arises 
from the difficulty of computing the matrix elements ( H2l) of the three-body Coulomb Green's 
function operator with sufficient numerical accuracy. Actually, the integrand includes os- 
cillatory functions whose amplitude grows very rapidly as the indices of the basis func- 
tions increase. As an example, Figure [2] shows the matrix element G^^q. ^, g (— |, with 
n = n' = 20 and S = ■y + i^e*'^, = — |. A comparison with Figure |3] shows that a relatively 
small change in the value of the angle can produce a large change in the amplitude. Note 
that in order to evaluate the matrix elements G^^q. q ( 1371) of the two-dimensional Green' 
function we resort to quadruple length arithmetic. 

The coefficients \ found to exhibit oscillations whose amplitude grows as in- 

creases. A simple way to damp the oscillations, is to multiply each matrix element [V]^,^ ^, 
by the Lanczos smoothing factors 

^ l-exp{-[ain-N)/Nf} 

= l-exp(-«2) ' (^^^ 

which attenuate [V]f,^ ^, with large indices rij, n'-. Thus, [V]fy^ in the expansion (!32l) are 

replaced by 

[Vko.'Il<<- (56) 

The optimal value for the parameter a in fl55l) is a ~ 3. The convergence behavior of 
the coefficient [a]g {rij = rrij = 0, j = 1, 2, 3) as the number is increased is presented in 
Figure 4. Note that such a simple remedy allows one to achieve convergence in the framework 
of the algebraic approach to two-particle scattering problems provided that the long-range 
part of the Hamiltonian is included into the "free" Green's function. The results for the 
first few coefficients [a]fy; are shown in Table [H In this calculation we obtained adequate 
convergence by including 15 basis functions (I3T!) for the coordinates j = 1, 2, 3. 



V. CONCLUSION 



The three-body continuum has been treated in terms of the generalized parabolic coordi- 
nates. A back-to-back electron emission from helium atom has been chosen as an example. 
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The Lippmann-Schwinger type equation for the continuum-state wave function has been 
solved numerically within the framework of the parabolic Sturmians approach. 

The potential in the basic integral equation is represented by the non-orthogonal part Di 
of the kinetic energy operator. Di treated in terms of the Cartesian coordinates ri3 and 
is a bounded operator. However, the change of variables ri3, — )■ ^j, rjj, j = 1, 2, 3 ([1]) 
transforms this operator into an unbounded one. Actually, the kernel of the operator Di 
grows without bound in the regions Qi and On the other hand, the triangle inequality 
for the vectors ri3, r23 is violated in these regions and therefore there does not exist a region 
in the real Cartesian coordinate system which corresponds to fii or ^2- Hence, this problem 
is cured by defining the kernel of Di in the regions fli and fl2 in an appropriate way (e. g., 
by setting the kernel equal to zero in these regions). For this purpose, we have approximated 
f 13 ■ ri2 and f 23 ■ ri2 which appear in the expression for Di by products of the projections of 
the vectors f^^ in the direction ki2. 

Besides, we have used the so-called outgoing approximation which assumes that the 
sought-for solution depends only on ^j, j = 1, 2, 3. Then the resulting integral equation 
compactness has been tested numerically by utilizing the parabolic Sturmian basis repre- 
sentation of the potential. In particular, convergence of the expansion coefficients of the 
solution has been obtained as the size of the basis is enlarged. 
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Figure 1: The path of integration C. The solid line is the part of C which remains on the physical 
sheet. The part of C which moves onto the unphysical sheet is shown by the dashed line. 
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Figure 2: The matrix element G,^ q. n' o ^) with n = n' = 20 and m = m' = on the 
contour C in Figure [1] for 99 = — ^ . 
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Figure 4: The convergence behavior of [ajg as the number N of basis functions (j3ip for each of 
the coordinates ^i, ^2, £,3 is increased. The coefficient values obtained without and with smoothing 
factors are denoted by squares and triangles, respectively. 
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Table I: Convergence of the first few "reduced" coefficients [a\syi- N is the number of the basis 
functions ^ (N = Ni = N2 = iVs) in the expansion ([32]) . 



N 




ni = n2 = 


= 713 = 




ni = 1, 7i2 




n3 = 




ni = 7i2 = 


= 0, 713 = 1 





1.00000000 + fO.OOOOOOOO 


0.33332667 - 


il.33332667 


-1.33333333 + iO.33333333 


(J 


1 
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iO. 16338707 





06484035 - 


il 


03904748 


-0 
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93751779 - 


iO. 00481271 





21244616 - 


il 


16523626 


-1 


09329190 


+ i0.31413236 


7 





85258743 - 


iO. 04143155 





28687445 - 


il 
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-1 


09100360 
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8 





88418841 - 


iO. 08046765 





25361563 - 


il 


06748525 


-1 


04328109 


+ i0.23545010 
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n 
u 




iU.uozy / oiD 


n 
u 


LLood t UO — 


6i 




1 

— 1 




-f- fcU.Z / DoOD ( t 


10 





90101510 - 


iO. 02409631 





23719582 - 


il 


12705918 


-1 


10091006 


+ i0.27238506 


11 





88184853 - 


iO. 02083012 





25897978 - 


il 


13142201 


-1 


11043846 


+ i0.25184571 


12 





87486976 - 


iO. 02519219 





26812774 - 


il 


12658516 


-1 


10795951 


+ iO.24433839 


13 





87207818 - 


i0.02578117 





27172914 - 


il 


12525415 


-1 


11016454 


+ i0.24216317 


14 
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iO. 02592191 





27617970 - 


il 
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-1 


11300765 


+ i0.23690466 


15 
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f0.02795984 
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il 
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-1 
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+ i0.23154456 


16 





86200713 - 
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